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SYNCHRONISATION AND ESTIMATION OF A TRANSMISSION CHANNEL 

The present invention relates generally to the field of digital communications; 
and, more particularly, to an efficient channel estimation method and 
arrangement for a digital communications system, such as a GSM/EDGE 
digital communications system, which provides a joint optimization of 
synchronization and equalizer window sizing. 

In a digital communications system, such as TDMA1, GSM, EDGE, or the 
like, data is encoded into symbols, packed into bursts and modulated prior to 
transmission via a physical transmission channel. At the receiver, the 
demodulation is performed by an equalizer which compensates for channel 
distortions, e.g. due to Inter-symbol interference (ISI). 

Inter-symbol interference (ISI) is an important problem in digital 
communications systems, including those systems which operate in 
accordance with the Global System for Mobile Communications (GSM) 
specifications. ISI is caused by time dispersion in the transmission channel 
over which a signal is transmitted, and adversely affects the quality of the 
received signal. In effect, ISI causes distortion of the transmitted signal 
which, in turn, causes symbol errors in the received signal such that it 
becomes difficult for the receiver to determine what data was actually sent. 

As is well-known in the art, the usual way to compensate for ISI in a/3SM 
system is to provide a channel estimation-based equalizer in the receiver. 
Basically, a model or estimate of the propagation channel over which a 
received signal was transmitted is created; and the equalizer then uses that 
information to estimate the sending symbols that best fit the received signal. 



30 



Correct burst synchronization and channel estimation, including estimation of 
the channel delay spread, are essential for the performance of the equalizer 
and other baseband processing. Burst synchronization and channel 
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estimation utilize the knowledge of a training sequence of symbols which is 
sent as part of the bursts, e.g. as a mid-ample in a burst. 

EDGE (Enhanced Data rates for Global Evolution) is an interface mode 
5 which has recently been developed for GSM Networks. EDGE'S principal 
features include new modulation and coding schemes which increase data 
capacity and speed in the air interface. EDGE is fully based on GSM and 
uses the same TDMA (Time Division Multiple Access) frame structure as 
GSM, such that it allows GSM operators to use existing GSM radio bands to 
1 0 offer wireless multimedia-based services and applications. 

In GSM/EDGE systems, the performance of the equalizer in combating ISI 
depends heavily on the quality of the channel estimation; and the quality of 
the channel estimation is, in turn, highly sensitive to the accuracy of burst 
15 synchronization (the term "synchronization" in this document signifies ,? burst 
synchronization") and the size of the equalizer window. 

To effectively combat ISI, the size of the equalizer window must be large 
enough to cover the maximum delay spread of the channel. However, an 

20 over-sized equalizer window will result in performance losses in channels 
with shorter delay spread due to inaccurate channel estimation, inadequate 
synchronization and increased noise contribution. The size of the equalizer 
window may be expressed as the number of filter taps of a corresponding 
equalizer filter. In the following, the size of the equalizer window will also be 

25 referred to as equalizer span or window span of the equalizer. In addition, 
synchronization can be optimized only with the knowledge of the equalizer 
window to capture maximum energy in the window span. In general, the 
interdependency of synchronization and equalizer window sizing makes 
efficient optimization difficult to achieve under different channel conditions. 

30 Joint optimization was deemed as too expensive in practical implementation. 

Because of this interdependency, a careful compromise has to be made for a 
traditional equalizer with fixed window span so that the size of the window 
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can provide adequate performance in long delay spread channels, such as 
Hilly Terrain (HT), without losing too much in short delay spread channels, 
such as Typical Urban (TU). 

5 When the equalizer window size is set, there are two approaches for burst 
synchronization, a correlation-based approach and a least square error (LSE) 
approach. As will be described more fully hereinafter, in current GSM 
receivers, synchronization is done by a correlation-based approach in which 
the central 16 symbols of a known 26-symbol training sequence is correlated 
10 with the central 16 symbols of the training sequence in the received signal. 

The synchronization is determined by either the center of gravity of the 
correlation or the maximum correlation energy in the equalizer window. In the 
LSE approach, for each possible synchronization point, an estimation of the 
15 channel is made, and synchronization is determined by comparing the 
expected and received training sequence via a least square error criterion. 

Neither of these approaches, however, is fully satisfactory. The correlation- 
based algorithm suffers from performance degradation due to inaccurate 
20 synchronization, especially in long dispersive channels, while a 
straightforward LSE-based algorithm suffers from a high degree of 
computational complexity (mainly due to multiple channel estimation). 

In particular, when non-white interferences, e.g. co-channel and adjacent- 
25 channel interferences as in GSM/EDGE systems, exist or when device 
related distortions exist, such as a DC off-set in homodyne receivers, 
correlation based algorithms suffer significant performance degradation. 

The above and other problems are solved by a method of estimating a 
30 transmission channel in a digital communications system including a receiver 
for receiving signal bursts of a communications signal via the transmission 
channel, the receiver comprising a channel estimation-based equalizer 
having an equalizer window, the method comprising 
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receiving a signal burst of the communications signal over the transmission 
channel, the signal burst including a sequence of training symbols; 

5 determining a desired synchronization position of the sequence of training 
symbols with respect to the received signal burst. 

The method is characterized in that the step of determining a desired 
synchronization position further comprises 

10 

determining a number of estimates of the transmission channel as a function 
of the synchronization position and a size of the equalizer window; 

determining a desired synchronization position and a desired size of the 
15 equalizer window by calculating an error measure based on the received 
signal burst and the determined estimates for a number of selected values of 
the synchronization position and of the size of the equalizer window, where 
the values of the size of the equalizer window are selected between 
predetermined upper and lower bounds. 

20 

Hence, as the desired synchronization position and the desired size of the 
equalizer window are determined based on a comparison of the received 
signal burst with a number of determined estimates for different values of the 
synchronization position and the size of the equalizer window, a joint 
25 optimization of synchronization and equalizer span is achieved and, thus, an 
improved performance of the equalizer for numerous channel conditions. 

By selecting the values of the size of the equalizer window between 
predetermined upper and lower bounds an interval for the span adaptation is 
30 defined, thereby limiting the necessary calculation steps. In the following, this 
interval will also be referred to as an aperture for the span adaptation. 
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Preferably, the predetermined upper and lower bounds define an interval 
having a wi'dth which is less than four channel taps, preferably less than 
three channel taps. Hence, the optimization of the equalizer span is limited to 
a narrow interval, thereby reducing the required computational resources. 

5 

It has been realized that the physical channel in a digital communications 
system does not change instantly with respect to the delay spread in contrast 
to, for example, the strength and phase of a signal which may experience 
rapid fluctuations. Based on this understanding, the aperture of the equalizer 
10 span can be made adaptive, i.e. the aperture for a current span optimization 
may be determined based on the determined equalizer spans for one or 
more previous bursts. 

In a preferred embodiment of the invention, the method further comprises 
15 determining the upper and lower bounds based on at least a desired size of 

the equalizer window as determined for a previously received signal burst. 

Consequently, the values of the size of the equalizer window are selected in 

a floating aperture with adaptively determined upper and lower bounds. 

Hence, the aperture for the span adaptation is made adaptive, and the 
20 predetermined upper and lower bounds are determined adaptively. 

Preferably, the adaptive determination of the floating aperture for the 
equalizer window further comprises filtering the values of the equalizer 
window sizes of previously received signal bursts. 

25 

In a further preferred embodiment, the step of determining the upper and 
lower bounds further comprises 

- determining a reference window size by means of an auto regressive filter 
which is updated for at least a subset of received bursts based on the 

30 determined desired size of the equalizer window for each of at least the 
subset of received bursts; and 

- selecting the upper and lower bounds to define an interval having a 
predetermined width around the determined reference window size. 
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Hence, the interval position is determined by a filtered version of the 
determined equalizer span, e.g. at predetermined time intervals such as for 
every received burst or at least for a subset of received bursts. The width of 
5 the interval may be limited to define a narrow aperture, e.g. 2-3 channel taps 
wide, in order to further reduce the computational effort. 

It is an advantage of the invention that it provides a joint optimization of 
synchronization and equalizer span that only requires little computational 
10 resources. 

It is a further advantage of the invention that it suppresses abrupt bursts-to- 
burst changes in the determined size of the equalizer window. Consequently, 
the effect of so-called "bad samples" in the channel estimation based 
15 approach is effectively reduced, thereby enhancing the robustness of the 
joint optimization against colored interferences and device distortions. 

The error measure may be a suitable measure comparing the received signal 
width an estimated signal as a function of the window size and the 
20 synchronization position. For example, the error measure may comprise the 
estimated noise power calculated based on the received signal and the 
estimated signal. Correspondingly, the step of determining the desired 
synchronization position and the desired size of the equalizer window size 
may comprise a Least Square Error optimization of the error measure. 

25 

In a preferred embodiment of the invention, the step of determining a desired 
synchronization position and a desired size of the equalizer window 
comprises calculating a difference between a calculated signal power of the 
received signal and a scaled power of a set of channel taps of the estimate of 
30 the transmission channel. 
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Consequently, an efficient calculation of the error function is provided which 
allows a calculation of the noise power without the need for calculating the 
noise signal. 

5 In yet another preferred embodiment of the invention, the method further 
comprises scaling the calculated difference with a penalty factor penalizing 
large equalizer windows. 

Preferably, the penalty factor comprises an exponential function of the 
1 0 equalizer span with a base determined by a length of the training sequence. 

Hence, it is an advantage of the invention that it reduces the risk of overfitting 
to noise when optimizing the equalizer span and synchronization position. 

15 Preferably, the communications signal comprises a signal in accordance with 
the GSM specifications or in accordance with the EDGE specifications. 

The present invention provides a more efficient channel estimation method 
by recognizing and utilizing two previously unexploited properties of 

20 GSM/EDGE training sequences which arise from the cyclic prefix structure of 
the training sequences. In particular, the present invention recognizes that 
the 26-symbol GSM/EDGE training sequences are, within certain ranges, 
both shift invariant and order invariant. The property of being shift invariant 
enables a channel estimation to be carried out with delayed (shifted) training 

25 sequence segments. This, in turn, permits: I. The ISI corrupted leading 
training sequence symbols to be avoided as much as possible in Ipng 
dispersive channels; and II. The leading taps of a channel to be estimated 
using the same training sequence segments, regardless of the size of the 
equalizer window. The property of being order invariant permits channels 

30 with different time dispersion, from 1-8 symbols (as complex polynomials of 
1-8th order) to be estimated without matrix inversion if any consecutive 16- 
symbol segment of the training sequence is used. 
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Hence, in a preferred embodiment, the step of determining a number of 
estimates of the transmission channel comprises selecting a subset of the 
number of training symbols. This results in a significant reduction in 
computational complexity in channel estimation for varying window size 
5 which can be adaptive to the actual delay spread of the channel. 

Furthermore, combining these two properties, a two-dimensional recursive 
relation can be established for different synchronization points and tap 
positions. 

10 

In a preferred embodiment, the step of determining a number of estimates of 
the transmission channel comprises determining channel taps of the estimate 
of the transmission channel for different synchronization positions and 
different sizes of the equalizer window by a two-dimensional recursive 
15 computation step. 

Consequently, the present invention provides a practical LSE joint 
optimization of synchronization and equalizer window sizing. 

20 It is a problem of the above prior art method that in channels with strong 
interferences a performance loss due to synchronization errors is more 
significant for the equalizer in a mobile receiver. Due to the limited length of 
the training sequences e.g. in GSM system, reliable synchronization is 
difficult to obtain. 

25 

Hence, according to yet another preferred embodiment of the invention, the 
step of determining a desired synchronization position and a desired size of 
the equalizer window by calculating an error measure based on the received 
signal burst and the determined estimates for a number of selected values of 
30 the synchronization position and of the size of the equalizer window further 
comprises selecting the values of the size of the synchronization position 
between predetermined upper and lower bounds; and the method further 
comprises determining the upper and lower bounds for selecting the values 
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of the size of the synchronization position based on at least a desired 
synchronization position determined for a previously received signal burst. 

Consequently, by using a narrow synchronization window centred around a 
5 reference synchronization position, synchronization is improved significantly, 
since fewer error-prone trials are necessary. The location of this narrow 
window is set with the information from synchronization positions of a 
previous burst, and thus "floats" depending on, e.g. the physical location of 
the mobile terminal. Hence, the candidate values of synchronization positions 
10 are selected in a floating aperture of adaptively determined upper and lower 
bounds. 

In a further preferred embodiment of the invention, the step of determining 
the upper and lower bounds for selecting the values of the size of the 
1 5 synchronization position further comprises 

- determining a reference synchronization position by means of an auto 
regressive filter which is updated for at least a subset of received bursts 
based on the determined desired synchronization position for each of at 
least the subset of received bursts; and 
20 - selecting the upper and lower bounds for selecting the values of the size 
of the synchronization position to define an interval having a 
predetermined width around the determined reference synchronization 
position. 

25 Hence, the adaptive determination of the floating aperture for the burst 
synchronization position further comprises filtering the values of the 
synchronization position of previously received signal bursts. Consequently, 
a simple and robust method of providing a floating reference position is 
provided. 

30 

In general, the present invention provides an efficient channel estimation for 
GSM/EDGE systems which, beside direct utilization, is fundamentally 
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important to achieving a practical implementation of joint LSE 
synchronization and equalizer window sizing. 

The present invention can be implemented in different ways including the 
5 method described above and in the following, an arrangement, and further 
methods and product means, each yielding one or more of the benefits and 
advantages described in connection with the first-mentioned method, and 
each having one or more preferred embodiments corresponding to the 
preferred embodiments described in connection with the first-mentioned 
10 method and disclosed in the dependant claims. 

It is noted that the features of the method described above and in the 
following may be implemented in software and carried out in a data 
processing system or other processing means caused by the execution of 
15 computer-executable instructions. The instructions may be program code 
means loaded in a memory, such as a RAM, from a storage medium or from 
another computer via a computer network. Alternatively, the described 
features may be implemented by hardwired circuitry instead of software or in 
combination with software. 

20 

The invention further relates to an arrangement for estimating a transmission 
channel in a digital communications system including a receiver for receiving 
signal bursts of a communications signal via the transmission channel, the 
receiver comprising a channel estimation-based equalizer having an 
25 equalizer window, the arrangement comprising 

means for receiving a signal burst of the communications signal over the 
transmission channel, the signal burst including a sequence of training 
symbols; 

30 

means for determining a desired synchronization position of the sequence of 
training symbols with respect to the received signal burst; 
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characterized in that the means for determining a desired synchronization 
position is adapted to 

determine a number of estimates of the transmission channel as a function of 
5 the synchronization position and a size of the equalizer window; and to 

determine a desired synchronization position and a desired size of the 
equalizer window by calculating an error measure based on the received 
signal burst and the determined estimates for a number of selected values of 
10 the synchronization position and of the size of the equalizer window, where 
the values of the size of the equalizer window are selected between 
predetermined upper and lower bounds. 

Hence, preferably, the arrangement further comprises means for determining 
1 5 said predetermined upper and lower bounds. 

The term means for receiving a signal burst includes any device or circuitry 
suitable for receiving signal bursts of the communications scheme used in 
the digital communications system. 

20 

The terms means for determining the desired synchronization position and 
means for determining said predetermined upper and lower bounds comprise 
any suitable processing means, where the term processing means includes a 
microprocessor, an application-specific integrated circuit, or another 
25 integrated circuit, a smart card, a general purpose computer adapted by 
suitable software, or the like. 

The invention further relates to a receiver for receiving signal bursts of a 
communications signal via a transmission channel in a digital 
30 communications system, the receiver comprising a channel estimation-based 
equalizer having an equalizer window and an arrangement for estimating the 
transmission channel, the arrangement comprising 
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means for receiving a signal burst of the communications signal over the 
transmission channel, the signal burst including a sequence of training 
symbols; 

5 means for determining a desired synchronization position of the sequence of 
training symbols with respect to the received signal burst; 

characterized in that 

1 0 the means for determining a desired synchronization position is adapted to 

determine a number of estimates of the transmission channel as a function of 
the synchronization position and a size of the equalizer window; and to 

15 determine a desired synchronization position and a desired size of the 
equalizer window by calculating an error measure based on the received 
signal burst and the determined estimates for a number of selected values of 
the synchronization position and of the size of the equalizer window, where 
the values of the size of the equalizer window are selected between 

20 predetermined upper and lower bounds. 

A receiver according to the invention may be part of an electronic equipment 
where the term electronic equipment includes all stationary and portable 
radio communication equipment and other handheld or portable devices. The 
25 term portable radio communication equipment includes all equipment such as 
mobile telephones, pagers, communicators, i.e. electronic organisers, smart 
phones, personal digital assistants (PDAs), handheld computers, or the like. 

It is a problem of the above prior art method that in channels with strong 
30 interferences a performance loss due to synchronization errors is more 
significant for the equalizer in a mobile receiver. Due to the limited length of 
the training sequences e.g. in GSM system, reliable synchronization is 
difficult to obtain. 
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Hence, according to another aspect of the invention, the above and other 
problems are solved when a method of estimating a transmission channel in 
a digital communications system including a receiver for receiving signal 
5 bursts of a communications signal via the transmission channel, the receiver 
comprising a channel estimation-based equalizer having an equalizer 
window, the method comprising 

receiving a signal burst of the communications signal over the transmission 
1 0 channel, the signal burst including a sequence of training symbols; 

determining a desired synchronization position of the sequence of training 
symbols with respect to the received signal burst from a number of estimates 
of the transmission channel corresponding to a number of selected values of 
15 the synchronization position 

is characterized in that 

the step of determining a desired synchronization position further comprises 
20 determining an upper and a lower bound for the synchronization position 
based on at least a desired synchronization position determined for a 
previously received signal burst; and selecting the number of selected values 
of the synchronization position between the determined upper and lower 
bounds for the synchronization position. 

25 

Consequently, by using a narrow synchronization window centred around a 
reference synchronization position, synchronization is improved significantly, 
since fewer error-prone trials are necessary. The location of this narrow 
window is set with the information from synchronization positions of a 
30 previous burst, and thus "floats" depending on, e.g. the physical location of 
the mobile terminal. 



WO 2004/047389 



PCT/EP2003/012202 



14 

In a preferred embodiment of the invention, the step of determining an upper 
and a lower bound for the synchronization position further comprises 

- determining a reference synchronization position by means of an auto 
regressive filter which is updated for at least a subset of received bursts 

5 based on the determined desired synchronization position for each of at 
least the subset of received bursts; and 

- selecting the upper and lower bounds for selecting the values of the size 
of the synchronization position to define an interval having a 
predetermined width around the determined reference synchronization 

10 position. 

Hence, a simple and robust method of providing a floating reference position 
is provided. 

15 Further preferred embodiments and the corresponding advantages 
correspond to the preferred embodiments of the first-mentioned method and 
their corresponding advantages. 

The invention further relates to an arrangement for estimating a transmission 
20 channel in a digital communications system including a receiver for receiving 
signal bursts of a communications signal via the transmission channel, the 
receiver comprising a channel estimation-based equalizer having an 
equalizer window, the arrangement comprising 

25 means for receiving a signal burst of the communications signal over the 
transmission channel, the signal burst including a sequence of training 
symbols; 

means for determining a desired synchronization position of the sequence of 
30 training symbols with respect to the received signal burst from a number of 
estimates of the transmission channel corresponding to a number of selected 
values of the synchronization position; 
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characterized in that 

the means for determining a desired synchronization position is further 
adapted to determine an upper and a lower bound for the synchronization 
5 position based on at least a desired synchronization position determined for a 
previously received signal burst; and to select the number of selected values 
of the synchronization position between the determined upper and lower 
bounds for the synchronization position. 

10 The invention further relates to a receiver for receiving signal bursts of a 
communications signal via a transmission channel in a digital 
communications system, the receiver comprising a channel estimation-based 
equalizer having an equalizer window and an arrangement for estimating the 
transmission channel, the arrangement comprising 

15 

means for receiving a signal burst of the communications signal over the 
transmission channel, the signal burst including a sequence of training 
symbols; 

20 means for determining a desired synchronization position of the sequence of 
training symbols with respect to the received signal burst from a number of 
estimates of the transmission channel corresponding to a number of selected 
values of the synchronization position; 

25 characterized in that 

the means for determining a desired synchronization position is further 
adapted to determine an upper and a lower bound for the synchronization 
position based on at least a desired synchronization position determined for a 
30 previously received signal burst; and to select the number of selected values 
of the synchronization position between the determined upper and lower 
bounds for the synchronization position. 
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The above and further aspects of the invention will be described in more 
detail in connection with a preferred embodiment and with reference to the 
drawing in which 

5 fig. 1 schematically illustrates the structure of a TDMA frame used in a digital 
communications system according to the GSM/EDGE standard; 

fig. 2 schematically illustrates the task of burst synchronization in a digital 
communications system; 

10 

fig. 3 is a schematic block diagram of a communications system according to 
an embodiment of the invention; 

fig. 4 is a flow diagram of burst synchronization and span adaptation 
1 5 according to an embodiment of the invention; 

fig. 5 illustrates the determination of the upper and lower bounds of a floating 
aperture for the equalizer span optimization; 

20 fig. 6 illustrates a two-dimensional recursive relation for the calculation of 
channel taps according to an embodiment of the invention; 

fig. 7 is a flow diagram of burst synchronization and span adaptation 
according to another embodiment of the invention; and 

25 

fig. 8 illustrates the determination of the upper and lower bounds of a floating 
aperture for the synchronization position. 

Fig. 1 schematically illustrates the structure of a TDMA frame 100 used in a 
30 mobile telecommunications system according to the GSM standard. Frame 
100 has a duration of 4.615 ms and accommodates 8 information channels 
(time slots) 101, designated 0-7 in Fig. 1. Each of the 8 time slots has a 
duration of 0.577 ms and contains a 148-bit signal portion and a guard 
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portion (not shown) which functions to maintain separation between signals 
in adjacent time slots. The 148-bit signal portion is generally referred to as a 
normal burst and comprises a first 3-bit tail bit section 102, a first 57-bit 
coded data section 103, a first 1-bit housekeeping bit section 104, a 26-bit 

5 trailing sequence section 1 05, a second 1 -bit housekeeping bit section 1 06, a 
second 57-bit coded data section 107 and a second 3-bit tail bit section 108. 
The 26-bits comprising the training sequence section 105 in a GSM digital 
communications system, is typically viewed as being divided into a centrally 
located portion 110 of 16-bits (sometimes referred to as the mid-amble 

10 portion), and side portions 109. and 111 of 5-bits each. Alternatively, in a 
GSM/EDGE digital communications system, the 26-bit training sequence 105 
can be also viewed as including a 16-bit white sequence and a 10-bit cyclic 
prefix. 

It is noted that, alternatively, other burst structures providing a suitable 
1 5 sequence of training symbols may be used. 

As mentioned previously, inter-symbol interference can occur in digital 
communications systems as a result of time dispersion in the transmission 
channel over which a signal is transmitted. In GSM/EDGE systems, the usual 
20 way to compensate for ISI is to provide a channel estimation based equalizer 
in the receiver. A correct burst synchronization is essential to the 
performance of the equalizer. 

Fig. 2 schematically illustrates the task of burst synchronization in a digital 
25 communications system. Fig. 2 shows a schematic example of a received 
burst 201 which is shown as the received power as a function of time. The 
shape of the burst is the result of the transmitted burst and the distortions 
introduced by the transmission channel, such as noise, multi-path 
propagation, etc. The task of burst synchronization involves the task of 
30 determining a suitable position p (203) of the training sequence 202 within 
the received burst. 
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As mentioned above, in known receivers, the equalizer window size m is set 
to a predetermined value and burst synchronization is performed by a 
correlation-based approach. By correlating the known training sequence t 
(202) with the received signal r (201), for all possible positions k within a 
5 synchronization window w, the following correlation coefficients may be 
calculated 



10 Here, the correlation length n is determined by the length of the training 
sequence 202. For GSM/EDGE, for example, the equalizer window size is 
often set to m=5 and the 16-symbol mid-ample in a 26-symbol training 
sequence is frequently used, i.e. n=21. With this approach, synchronization 
may be achieved by a maximum tap energy approach, i.e. the 

1 5 synchronization position p may be determined as 



20 Another approach utilizes a center of gravity condition, according to 



In an alternative approach, burst synchronization may be achieved by a 
Least-Squares Error (LSE) method. This approach requires an initial channel 
estimation h p for each possible synchronization position p. The desired 
30 synchronization position is then determined by minimizing a square error 
according to 



c(k)=£r(i+k).t(i), k-0,1 w-1. 
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p = argmax £|c(k + i)| 
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(0) 



m-1 

fp(k)=£h p (m-i-i)t(k+i) 

1=0 

is an estimate of the received symbol based on the initial channel estimate 
hp. 

10 However, neither of these approaches is fully satisfactory. The correlation- 
based algorithm suffers from performance degradation due to inaccurate 
synchronization, especially in long dispersive channels, while a 
straightforward LSE-based algorithm suffers from a high degree of 
computational complexity (mainly due to multiple channel estimation). 

15 

In particular, when non-white interferences, e.g. co-channel and adjacent- 
channel interferences as in GSM/EDGE systems, exist or when device 
related distortions exist, such as a DC off-set in homodyne receivers, 
correlation based algorithms suffer significant performance degradation. 

20 

Further, when the channel delay spread is different from the assumed value, 
considerable performance losses may result, not only because an 
inadequate span degrades the equalizer performance, but also because an 
inadequate synchronization may occur. 

25 

Fig. 3 is a schematic block diagram of a communications system according 
to an embodiment of the invention. The communications system comprises a 
transmitter 301 and a receiver 303. A signal is transmitted from the 
transmitter 301 to the receiver 303 over a transmission channel 302. The 
30 transmission characteristics of the transmission channel frequently change 
due to a variety of factors, including the motion of the mobile terminals, the 
fluctuation of multi-path (time dispersive) propagation channels and variant 



where 

5 



WO 2004/047389 



PCT/EP2003/012202 



20 

interferences introducing noise. As is well-known to those skilled in the art, 
the multipath channel and the noise component can adversely affect the 
quality of the received signal, e.g., cause intersymbol interference; and this 
necessitates that the received signal be corrected, typically by means of 
5 channel estimation-based equalizer. 

Specifically, as shown in Fig. 3, at the receiver 303, the received signal, 
which includes the transmitted signal with intersymbol interference and noise, 
is, after being filtered by a filter 304, sent to an equalizer 308 which attempts 

10 to correct the received signal utilizing an estimate of the transmission 
channel. The equalizer 308 generates data representing an estimate of the 
actual transmitted symbols for further processing by the receiver. The 
receiver 303 further comprises a channel estimation nodule 307 for 
generating a channel estimate which is fed into the equalizer 308. The 

15 channel estimate, in turn, receives a burst synchronization position p and a 
window size m of the equalizer window as an input. 

According to the invention, the receiver 303 comprises a joint synchronization 
and span adaptation module 305 which performs a joint determination of a 
20 synchronization position p and a size m of the equalizer window to be fed into 
the channel estimation module 307. 

The joint optimization of the equalizer span and the synchronization position 
performs an optimization within a predetermined interval of window sizes. 

25 Accordingly, the receiver 303 further comprises an aperture module 306 
which determines an upper and a lower bound m u and mi, respectively, of the 
window size to be used in the optimization step for the subsequent burst. The 
aperture module 306 receives the determined window size m for a current 
burst from the span adaptation module 305 and feeds resulting bounds for 

30 the next burst back to the span adaptation module 305. In a preferred 
embodiment, the aperture module 306 further determines an upper and a 
lower bound p u and pi, respectively, of the synchronization position to be 
used in the optimization step for the subsequent burst. Hence, the aperture 
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module 306 receives the determined synchronization position for a current 
burst from the span adaptation module 305 and feeds resulting bounds for 
the next burst back to the span adaptation module 305. 

5 It is noted that in connection with an EDGE receiver, it is preferred to 
combine the receiver filter with a noise whitening filter in order to suppress 
interferences and, thus, to avoid the necessity of a narrow receiver filter. 

Fig. 4 is a flow diagram of burst synchronization and span adaptation 
10 according to an embodiment of the invention. Upon receipt of a signal burst 
r(k) in step 401, a Least Square Error algorithm is performed to determine the 
equalizer window span m and the synchronization position p as illustrated by 
the loop comprising steps 402-405. 

15 According to this embodiment, the equalizer window span m and the 
synchronization position p are determined as 

(P>nr0 opt =argmin(8 2 (p,m))=argmin(a m .e 2 (p,m)) P = °» 1 -- w " 1 (1) 
p.m p,m m-m|,... l m u 

20 where e(p,m)=a m/2 'e(p,m) is a generalized error measure derived from the 
error function 

e (P. m ) = ^Shk)-fp,m(kf- 



m-1 

W( k )=Z h P,m( m " 1 - i )' t ( k+1 ) 
i=0 



Here 

25 
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is an estimate of the received signal expressed as a function of m and p, i.e. 
e(p,m) is a measure of the noise power introduced by the estimated 
transmission channel h p , m on the training sequence t. 

5 If the noise power is used directly as an error function for optimization, the 
channel window size m tends to be larger than necessary. This is because 
the additional parameters of a larger filter could adjust themselves to 
particular features of the specific realization of the noise, i.e. the 
phenomenon sometimes referred to as overfitting. To avoid this, a penalty 

10 factor <x m is introduced according to the Minimum Description Length 
Principle (MDL) (see e.g. R. Johansson, "System Modeling and 
Identification", Prentice Hall, 1993). The penalty factor penalizes large 
equalizer spans m, and therefore suppresses the effect of overfitting. In a 
preferred embodiment, the factor a is determined by the length n of the 

15 training sequence ct= n 1/n . Hence, in the example of a 16-symbo! training 
sequence a=1.189. It is noted, however, that alternatively other functions 
penalizing large window spans may be used. 

The optimization of eqn. (1) is performed both for p and m, where m is varied 
20 in an interval between the upper and lower bounds m u and mi, respectively, 
as will be described in greater detail below. 

Still referring to fig. 4, the optimization loop thus comprises the following 
steps: 

25 

Step 402: Select a pair of values (p,m) within the intervals indicated in eqn. 
d). 

Step 403: Estimate the transmission channel h p , m for the selected values of p 
30 and m based on the stored training sequence t (408). An efficient method of 
determining the channel estimate for a training sequence of the GSM system 
will be described in connection with fig. 6. 
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Step 404: Calculate the generalized error e 2 (p,m). As will be shown in greater 
detail in connection with fig. 6, for 16-symbol GSM training sequences and 
with m<8, the generalized error e 2 (p,m) may be expressed as 

5 e 2 (p f m) = a m (|r| 2 -16|h Pim | 2 ), s 2 >0. 

Hence, in this case, the error measure e 2 (p,m) may be efficiently calculated 
as the difference of the received signal power |r| 2 and the scaled power of the 
estimated channel taps 16.|h p , m | 2 , i.e. without the need for calculating the 
10 actual noise. It is noted, that the constraint e 2 >0 is explicitly enforced, since 
the expression |r| 2 - 16|h p<m | 2 may become negative at incorrect 
synchronization positions. Hence, when 16|h p , m | 2 £ |r| 2 no further calculation 
is necessary and the hypothetical synchronization position is rejected 
immediately, since the constraint is not met. 

Step 405: Repeat the above steps until a suitable minimum is found. In one 
embodiment, the error is calculated for all possible pairs (p,m) within the 
above intervals, and the minimum of all calculated values is determined as 
the optimal set of values (p,m) op t. 

20 

When the optimal set of values (p,m) op t is determined, in step 406 the 
aperture for the span optimization of the subsequent burst, i.e. the upper and 
lower bounds m u and mi, is determined and stored in a memory 407 for use 
in the subsequent optimization. According to this embodiment, an adaptive 
25 aperture is achieved by use of a simple auto regressive (AR) filter with one 
state m s that is updated at each burst: 

m 8 (t) = a.m s (t-1) + b.m(t) l 

30 where m(t) is the current equalizer span determined by the above 
optimization. Hence, the state m s (t) corresponds to a weighted mean of the 
previously determined equalizer spans, where the spans determined for the 
0 most recent bursts are weighted strongest. A suitable initial value is 
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assumed, e.g. m s (t=1) = m(t=1) = (m m in + m ma x)/2. The relative weights are 
determined by the parameters a and b. Preferably a and b are selected in the 
interval [0,1], preferably such that a>b, e.g. aefO.8,0.9] and be[0.1,0.2], for 
example a= 0.875 and b=0.125. Larger values of a reduce the relative 
5 influence of the most recent changes and vice versa, i.e. preferably the 
parameters a and b are selected according to the time scale of the variation 
in the delay spread. 

The aperture of the equalizer span optimization is then determined as: 

10 

mi= round(m s (t-1)) - Wi, rri|£ m m i n 
m u = round(m s (t-1 )) + W 2 , m max ^ m u , 

i.e. mi and m u are lower and upper bounds, respectively, of an interval 
15 around m s . The constants Wi and W 2 determine the width of the aperture. 
For example, Wi=W 2 =1 yields an aperture that is 2-3 taps wide. 

The constant upper and lower limits m ma x and m m in ensure that mi and m u lie 
within desired limits, e.g. m m in = 4 corresponding to typical urban (TU), and 
20 m max = 8 corresponding to hilly terrain (HT). 

Fig. 5 illustrates the determination of the upper and lower bounds of a floating 
aperture for the equalizer span optimization. Lines 501 and 502 illustrate the 
values mi and m u , respectively, selected for different values of m s according 
25 to eqn. (3). Lines 503 and 504 illustrate an alternative constant choice of 
bounds. 

Fig. 6 illustrates a two-dimensional recursive relation for the calculation of 
channel taps according to an embodiment of the invention. The optimization 
30 method described in connection with fig. 4 above requires that w x (m u - mi + 
1) channel taps are estimated. In the following a computationally efficient 
channel estimation is described which exploits certain properties of the 
GSM/EDGE training sequences. 
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Hence, in the following a GSM/EDGE training sequence is considered: 

t = [toti ... Wi ...tn.i] H . 

5 

For real training symbols, the Hermitian conjugate [] H corresponds to simple 
matrix transposition. 

For n>m, a channel estimation can be obtained by a LSE algorithm according 
10 to 

H = O- 1 T H r, 

where r is the vector of (n-m) received symbols, delayed by m-1 symbols, 
15 and 



20 



25 



30 



0=T H T 

is the autocorrelation matrix of dimension m x m the m x (n-m) dimensional 
Toeplitz-like matrix T of the training sequence: 



T = 



*m-1 tm — *n-1 
tm-2 tm-1 - : 



L *0 U - t n - 



m. 



In GSM/EDGE systems, 8 training sequences are specified for a normal 
burst A transmission takes one of the sequences according to the so-called 
Training Sequence Code (see e.g. GSM 05.02, V8.0.0, ETSI, March 1999). 
Each training sequence consists of 26 symbols. 

The GSM/EDGE training sequences have two properties, not previously 
exploited, which permit channel estimation to be carried out more efficiently. 
More particularly, the autocorrelation matrix of any consecutive 16-symbol 
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segment of the GSM/EDGE training sequence is, within certain range, both 
shift invariant and order invariant. 

Property I. Shift Invariant: 

5 

<D(m,d) = T H (m,d) T(m,d)=1 6 l m , (26 -(16 + m)+1 ) ^ d ;> 0, 

where <£(m,d) and T(m,d) are written as functions of the number channel- 
taps m and the index d of the first training symbol in T, and where l m is the m- 
1 0 dimensional identity matrix. 

Property II. Order Invariant: 

<D (m,d)=T H (m,d) T(m,d) = 16 l m , 7 ^ m ^ 1. 

15 Property I is due to the cyclic prefix nature of the GSM/EDGE training 
sequences, in which the 26-bit training sequence is constructed by 
prepending the last 1 0 bits to its front. 

This property enables a channel estimation to be carried out on delayed 
20 (shifted) training sequence segments. This flexible delay brings two important 
benefits: 

(1) ISI corrupted leading symbols can be avoided in the computation; and 

25 (2) the leading taps of a channel can be estimated using the same training 
sequence segments, regardless of the size of the equalizer window. 

Property II enables channel estimation to be done for different equalizer 
windows ranging from 1 to 7 taps. Furthermore, it has been found that when 
30 m = 8, four of the eight specified training sequence code (TSC) satisfy the 
above order invariants, and the auto-correlation matrix for the other four 
training sequences is in a specific sparse form: 
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16 0 

0 16 



+ 4 



3>(8,d)= 



0 

±4 0 



16 



5 



16 



The inverse of <B(8,d) has the same sparse, symmetric structure: 



a 0 
0 1 



10 



a>-'(8,d)=-L 



0 

0 0 



1 



a 



15 where a = 1 .0667, B= ± 0.2667. This means that in the above estimation of 8- 
tap channels, only minor adjustments are necessary for certain TSCs: 



Property II indicates that not only a 5-tap channel, but all the 1-8 tap 
channels can be estimated efficiently without matrix inversion, if any 
consecutive 16-symbol segment of the training sequence is used. Utilization 
of these properties has a fundamental impact on the computational 
25 complexity of the channel estimation for different equalizer window sizes. 

When an LSE approach is used for joint synchronization and equalizer 
window sizing, for a hypothetical sync position p and an equalizer window 
size m, by utilizing the properties discussed above, channel estimation can 
30 be computed as 




20 



H = o- 1 T H r = (1/16)T H r. 
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10 



15 



20 



25 



With the exception of the constant factor, it can be further expressed in an 
explicit convolution form, 



This can be considered as a FIR-style computation. For the next tap with 
k+1 , a recurrent relation can be established: 

h P> m(k+1) = h p+ i, m (k) + t(m-1-k)-r(p+m)-t(m + (n-1)- k)-r(p+m + n). (5) 

Using this recursion in an IIR-style computation, a new tap can be calculated 
using only 4 real MAC operations (since the training sequences are real). 
This two-dimensional (sync-point tap-position) recursive relation can be 
depicted as shown in Fig. 6 for m = 5, w = 8. Computation of the last row and 
the last column (indicated by the area 601) is done by using equation (4), the 
rest can be recursively calculated using eqn. (5). 

To avoid channel estimation error caused by ISI corrupted receiving training 
symbols, by property I, a training sequence segment that is close to the end 
may be chosen. In practice, the last but one m + 15 symbols of the training 
sequence are chosen in the estimation (the final symbol in the training 
sequence is not used, to avoid possible precursor ISI corruption). 

As mentioned above, using the above properties, a simplified expression for 
the error measure e 2 may be achieved. To show this, the channel estimate 
hp, m and the corresponding prediction of the received signal may be written in 
matrix notation as follows: 



h p ,m (k) = J t(m + i - k) • r(p + m - i), 



0<p^ w 
0<k<m 



(4) 
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Thus, the error e 2 can be written as 

e 2 =(r-f) 2 = r H r - r H T(T H T)r 1 T H r. 
5 Using the above relation H = O" 1 T r = (1/16) T r, this can be simplified as 

e 2 = |r| 2 -16|h Pim | 2 , 

thereby yielding an efficient calculation of the squared error as a weighted 
10 energy difference of the 16-symbol received training sequence and the m-tap 
channel estimate, which can efficiently be calculated by the above described 
recursion. 

The presented method has been put into simulation and verified. Simulations 
15 show that a significant performance gain is observed in long time dispersive 
channels with joint optimization compared to the correlation based 
synchronization. The joint synchronization and equalizer span adaptation 
significantly improves the performance of the equalizer in numerous channel 
conditions, i.e. in the range of 1-4dB, with a comparable computational cost 
20 as the correlation based approach. 

As mentioned above, the adaptive aperture in equalizer span adaptation 
further enhances the robustness of the joint optimization against colored 
interferences and device distortions. Correspondingly, a further improvement 
25 of more than 1 dB has been observed in simulations. 

Fig. 7 is a flow diagram of burst synchronization and span adaptation 
according to another embodiment of the invention. The method comprises 
step 701 of receiving a signal burst r(k), and the method comprises the 
30 performance of a Least Square Error algorithm to determine the equalizer 
window span m and the synchronization position p as illustrated by the loop 
comprising steps 702-705. 
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According to this embodiment, the equalizer window span m and the 
synchronization position p are determined as 

(p,m) opt =argmin(8 2 (p l m))=argmin(a m .e 2 (p > m)) P = P| '"" p " ( 6 ) 
p,m p p m m-m|,...,m u 



where e(p,m)=a m/2 e(p,m) is a generalized error measure derived from the 
error function 



e (P. m ) = jZ| r ( k )- f P,m(k)| 2 . 

10 Here 



m-1 

fp,m( k )=S h P,m(m-1-i)^t(k + 1) 
i=0 



is an estimate of the received signal expressed as a function of m and p, i.e. 
15 e(p,m) is a measure of the noise power introduced by the estimated 
transmission channel h Pjm on the training sequence t. 

As described in connection with fig. 4, a penalty factor <x m is introduced that 
penalizes large equalizer spans m, and therefore suppresses the effect of 
20 overfitting. In a preferred embodiment, the factor a is determined by the 
length n of the training sequence a= n 1/n . Hence, in the example of a 16- 
symbol training sequence ct=1.189. It is noted, however, that alternatively 
other functions penalizing large window spans may be used. 



25 The optimization of eqn. (6) is performed both for p and m, where m is varied 
in an interval between the upper and lower bounds m u and mi, respectively, 
as was described in connection with figs. 4 and 5. Furthermore, according to 
this embodiment, p is varied in an interval between the upper and lower 
bounds p u and pi, respectively. 
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Still referring to fig. 7, the optimization loop thus comprises the following 
steps: 

5 Step 702: Select a pair of values (p,m) within the intervals indicated in eqn. 
(6). 

Step 703: Estimate the transmission channel h Pim for the selected values of p 
and m based on the stored training sequence t (708). An efficient method of 
10 determining the channel estimate for a training sequence of the GSM system 
has been described in connection with fig. 6. 

Step 704: Calculate the generalized error e 2 (p t m). It was described in 
connection with step 404 of fig. 4, how the error measure e 2 (p,m) may be 
1 5 efficiently calculated for 16-symbol GSM training sequences and with m<8. 

Step 705: Repeat the above steps until a suitable minimum is found. In one 
embodiment, the error is calculated for all possible pairs (p,m) within the 
above intervals, and the minimum of all calculated values is determined as 
20 the optimal set of values (p,m) 0 pt. 

When the optimal set of values (p,m) 0 pt is determined, in step 706 the 
aperture for the span optimization of the subsequent burst, i.e. the upper and 
lower bounds m u and mi, is determined and stored in a memory 707 for use 
25 in the subsequent optimization, as was described in greater detail in 
connection with step 406 of fig. 4. According to this embodiment, additionally 
an aperture for the optimization of the synchronization position of the 
subsequent burst, i.e. the upper and lower bounds p u and p (l is determined 
and stored (707) for use in the subsequent optimization. 

30 

Similarly to the determination of mi and m Ut an adaptive aperture for the 
synchronization position is achieved by use of a simple auto regressive (AR) 
filter with one state p s that is updated at each burst: 
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Ps(t) = A.p s (M) + B.p(t) I (7) 

where p(t) is the current synchronization position determined by the above 
5 optimization. Hence, the state Ps(t) corresponds to a weighted mean of the 
previously determined synchronization positions, where the positions 
determined for the most recent bursts are weighted strongest. A suitable 
initial value is assumed, e.g. p s (t=1) = p(t=1) = (p m in + Pmax)/2. The relative 
weights are determined by the parameters A and B. Preferably A and B are 

10 selected in the interval [0,1], preferably such that A>B, e.g. A<=[0.8,0.9] and 
Be[0.1,0.2], for example A= 0.875 and B=0.125. Larger values of A reduce 
the relative influence of the most recent changes and vice versa, i.e. 
preferably the parameters A and B are selected according to the time scale 
of the variation in the synchronization position. It is noted that the 

1 5 parameters A and B may be selected to be different from or equal to the 
corresponding parameters a and b, respectively, used in the determination of 
the equalizer span aperture. 

The aperture of the optimization of the synchronization position is then 
20 determined as: 

pi s round(p s (t-1 )) - w/2, (8) 
Pu= Pi + w, 

25 i.e. pi and p u are lower and upper bounds, respectively, of an interval of width 
w around p s . Hence, the synchronization position only needs to be 
determined within a narrow interval, e.g. having a width w of 3-5 symbols. It 
is noted that, in an alternative embodiment, a asymmetric interval around p s 
may be used instead. Furthermore, in one embodiment, the upper and lower 

30 bounds p u and pi may be restricted to lie between predetermined limits p mln 
and p max . 
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Hence, a simple auto-regressive filter is used to smooth the optimal position 
of the narrow window. 

It is noted that, in an alternative embodiment, the floating aperture of the 
5 synchronization position may be applied in connection with a synchronization 
according to eqn. (0>. 

A floating aperture for synchronization window not only improves equalizer 
performance, but also reduces the computational complexity since, at each 
1 0 burst, fewer synchronization positions need to be tested. 

Fig. 8 illustrates the determination of the upper and lower bounds of a floating 
aperture for the synchronization position. Lines 801 and 802 illustrate the 
values pi and p Ul respectively, selected for bursts at different times t 
15 according to eqns. (7) and (8). Lines 803 and 804 illustrate an alternative 
constant choice of bounds. 

The invention has primarily been described in connection with GSM/EDGE 
training symbols. However, the invention is not limited to GSM/EDGE, but 
20 may be applied to signal bursts of other communications schemes which 
utilize suitable training sequences for burst synchronization. 



